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Abstract: Inspired by the recently discovered non-relativistic fermionic fixed points, 
we investigate how the presence of bulk dipole coupling modifies the spectral function 
at one of these novel fixed points. As a result, although the infinite flat band is always 
visible in the presence of the bulk dipole coupling as well as chemical potential, the 
band is modified in a remarkable way at small momenta up to the order of magnitude 
of bulk dipole coupling. On the other hand, like a phoenix, a new Fermi surface sprouts 
from the formed gap when the bulk dipole coupling is pushed up further such as to 
overshadow the charge parameter, which is obviously different from what is found at 
the relativistic fixed points. 
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1. Introduction 

Arising from the string theory, AdS/CFT correspondence has provided us with a new 
paradigm to study strongly coupled many body phenomena by relating them to a 
single or few body problem in the bulk classical gravitational backgrounds with one 
extra dimension. In applying AdS/CFT correspondence, there are two approaches 
which have been taken. One is the so called bottom-up approach, where the bulk setup 
is devised in a simple way but the microscopic understanding of dual field theory is 
generically lacking. The other is the so called top-down approach, where the whole 
setup can be embedded in the more sophisticated string theory or M theory such 
that the microscopic content of boundary field theory is well understood, although 
the system may not be realized in Nature. However, such a shortcoming is somehow 
mitigated by focusing judiciously on the IR physics, where it is believed that the same 
universal behaviors can be extracted from this sort of holographic duality for those 
realistic systems which may have different microscopic contents. 

With this in mind, attempts to apply this holographic technique, with some suc- 
cess, have gone beyond high energy physics, especially to condensed matter physics in 
the recent years. Obviously, it is kind of win-win game to apply the holography to con- 
densed matter physics. On the one hand, there exist many strong coupled systems in 
condensed matter physics, which are intractable by the conventional approaches. While 
AdS/CFT correspondence, as kind of strong/weak duality, can offer unprecedented in- 
sight into the dynamics of these systems. On the other hand, unlike other disciplines 
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such as high energy physics and cosmology, only with some tabletop experiments, does 
condensed matter physics allow one to cook up matter such that various vacuum states 
and phases can be created in the laboratory. With this available landscape of man- 
made multiverse to scan, it is highly possible for us to hit some universality classes 
extracted from the holography. In this sense, condensed matter physics may in turn 
provide us with the first experimental evidence for AdS/CFT correspondence. For a 
review of this exciting subject, please refer to [0, |, H, |, ^. 

In particular, partly triggered by the lack of a proper field theoretical framework to 
explain for the mysterious behaviors of existing non- Fermi liquids, endeavors have been 
made recently to study the spectral function of fermionic operators by holographically 
manipulating the bulk Dirac field which is minimally coupled to gravity and gauge 
fields, where the emergence of Fermi surface has been identified with a rich spectrum 
of non- Fermi liquid behaviors[|], ^ H, H |lO[ 0. Later on, inspired by the generic top- 
down models, the bottom-up dipole coupling has been added to the previous minimally 
coupled Dirac field and its effects on the fermionic correlator have been investigated, 
where a possible dynamical gap opens up when the dipole coupling is tuned to be 
large enough |T^, |13], Generically, in the context of AdS/CFT correspondence, not 
only the modification of bulk dynamics, but also the change of boundary conditions 
can alter the boundary field theory. Actually, in the latter manner, the holographic 
non-relativistic fermionic fixed points have been implemented most recently by adding 
a Lorentz violating boundary term instead of the conventional Lorentz covariant one to 
the bulk minimally coupled Dirac action, where a holographic flat band is achieved [|17||. 

Along this line, this paper is intended to insist on one of these novel non-relativistic 
fixed points and investigate how the fermionic correlator is modified when we turn on 
the bulk dipole coupling. In the next section, we shall provide a brief review of how 
both of the relativistic and non-relativistic fermionic fixed points are implemented by 
holography. In Section ^, after building up the holographic framework to extract the 
fermionic correlator at the non-relativistic fermionic fixed point, we shall present our 
numerical results for the relevant quantities associated with the fermionic correlator in 
the presence of bulk dipole coupling. Conclusions and discussions will be addressed in 
the end. 

2. Holographic implementation of various fermionic fixed points 

Start with the bulk action for a probe Dirac fermion with the mass m, charge q and 

^It is noteworthy that the impact of bulk dipole coupling on fermionic correlator has also been fully 
explored from the purely top-down construction in and JlGt - 
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magnetic dipole coupling p 

Jm '-^ 
in the following fixed background, i.e., 

ds^ = r^[-f(r){dtf + {dx'f + [dx^f] + 1^^, A, = A^idt),. 



(2.1) 



(2.2) 



Here i, = ^tp*, '$ = {e^fT^'lda + \[{u:p<r)aVP'' - ig^J, and f = \V^''^{e^Y{e,fFa^,, 
where (e^)" form a set of orthogonal normal vector bases, and Gamma matrices satisfy 
I^At ^ 2?7^^ with the spin connection {uj^,y)a = (e^)bVa(e^)^ T^'' = l[T^',T•'], and 
the field strength F = dA. In addition, the emblackening factor and gauge potential 
are given by 



/ 



1 + , Q 



+ ^,At = gpQil- 



(2.3) 



r 

which arises as the charged black hole solution to the equation of motion following from 
the bulk action for the gauge field coupled to AdS gravity with the gauge coupling gp, 
i.e., 

1 
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M 
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(2.4) 



where L is the curvature radius, and has been set to unity along with the horizon 
radius. By holography, such a charged black hole places the probe fermion into the 
dual strongly coupled soup with the finite temperature and chemical potential as 



3-g2 

r= ,1^^ gpqQ- 



(2.5) 



Now to have a well-defined variational principle for the Dirac action, a boundary 
term must be added. To see this, let us make the variation of bulk action, which gives 
rise to 



^Sbulk = i / d'^xy/^ 



M 



'^S'iIj($ -m- ipf)'^ - i/j^^ + m + ipf)5'^ 



+- 



d^xV—h{ip-S'i/j+ — i/j+Si/j^ + — 5'0_'0+) 



(2.6) 



dM 



where h = — is the determinant of induced metric on the boundary, and ip± — 

Qrr 



^{1 ±T'^)ip. Note that the bulk Dirac equation 

0-711- ip^)ijj = 



(2.7) 
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is first order. So it is not allowable to fix all components of ^jj, namely both of ^p+ 
and ip-. Instead we must somehow fix simply half of the components of ip, which can 
actually be implemented by adding a boundary term. Moreover, it turns out that the 
variational principle can be achieved by adding a number of different boundary terms, 
depending on the specific value of the mass parameter m^. 

The conventional boundary term chosen to be added is the Lorentz covariant one, 

i.e., 

Sbdy = ^ / v^^^ [ V^(^-^+ + (2.8) 

^ JdM ^ JdM 

whereby the variation of the full on-shell action is given by 

SSbuik + SSbdy = I / d^xv^(V^-<5V'+ + 5^+^-), (2.9) 

JdM 

which indeed vanishes if and only if the Dirichlet boundary condition is imposed on 
This sort of choice of boundary condition is usually referred to as the standard 
quantization for the Dirac field. As a result, the dual boundary field theory is a Lorentz 
covariant CFT where loosely speaking plays a role of fermionic source, and the dual 
operator is given by with dimension 

A[V^_] = ^ + m. (2.10) 

When the mass parameter is tuned into the window < m < |, we can have other 
boundary conditions to choose. The first one is the so called alternative quantization, 
which can be implemented simply by adding the boundary term with opposite sign, 
i.e., 

Smv = --1 V^i^^ = -'- I ^{^-^+ + (2.11) 

^ JdM ^ JdM 

This now results in a well defined variational principle if and only if one imposes 
the Dirichlet boundary condition on The dual boundary field theory is still a 
Lorentz covariant CFT with the fermionic source and dual operator interchanged. The 
dimension of operator is thus given by 

A[^+]- = ^-m. (2.12) 

We have other boundary conditions to choose if we are not only interested in the Lorentz 
covariant boundary field theory. In particular, as shown in [rf|, a non-relativistic 



^For simplicity but without loss of generality, we shall focus on the case of m > in what follows. 
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fermionic fixed point can be implemented by adding the following Lorentz violating 
boundary term, i.e., 



Si 



1 



bdy 



(2.13) 



2 JdM 

To see explicitly what both of the corresponding fermionic source and dual operator 
look like, we would like to firstly choose our Gamma matrices once and for all as follows 

-(T^ 







-cr 



ia^ 
ia^ 



-a^ 






la^ 



(2.14) 



where a* are Pauli matrices^. With this kind of choice of Gamma matrices, we can 
express ip^ and ip^ as 







-h) 





w 



(2.15) 



Hereby the variation of bulk on-shell action can be massaged as 
SSbuik = -7^ / d^x{6zlyi+6zly2-6ylzi-6ylz2-zl6yi-zl6y2+yl6zi+yl6z2). (2.16) 

^ JdM 

On the other hand, hj ip = 4'+ + "^-5 the Lorentz violating boundary term can be 
casted into 

(2.17) 



1 



Sbdy = -7; I d-^x^zlyi + ylzi + z[y2 + y\z2). 

^ JdM 

Whence the variation of full on-shell action is given by 

^Siyuik + ^Sbdy = -\ I d^x[5{z\ + z\){yi + 2/2) + 5{y{ - yl){z2 - z\) 

^ JdM 



,t. 



+(4 - zl)6iyi - y2) + {y\ + yl)5{zi + Z2)] 
- [ d^x{5Z\Yi + Zl5Y2 + Yl5Zi + 5Y^Z2) 

JdM 



(2.18) 



where we have defined (1^1,1^2) = ^(1/1 + 2/2,1/1-1/2), and (^1,^2) = ■^{zi + Z2, Z2- Zi). 
Hence the corresponding fermionic source and dual operator are somehow given by 
(Zl, Y2) and (Yi, Z2) respectively, and the dimension of operator reads 



A[ri] = ^ + m,A[Z2] = ^ 



m. 



(2.19) 



•^Note that our choice of Gamma matrices is different from that in [|T^. We are believed that our 
choice is more suitable for performing the practical calculation, which can be seen later on. 
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Regarding other possible choices of boundary conditions and how they are related 
to one another through the Wilsonian RG flow by the double trace deformation, please 
refer to 



17, 1 



In the subsequent section, we will focus only on the above non- 
relativistic fixed point, to extract the boundary fermionic correlator from the bulk 
dynamics and see how the magnetic dipole coupling affects the fermionic correlator by 
numerical calculation. 



3. Holographic non-relativistic fermion with bulk dipole cou- 
pling 

3.1 Holographic setup 

Generically the relevant information regarding the fermionic system can be read out 
of its single particle fermionic correlator, namely the retarded Green function Gr. 
For example, the spectral function A{u, k), which can be measured experimentally by 
Angle Resolved Photoemission Spectroscopy(ARPES), is given by the imaginary part 
of TtGr. Now we shall show how such a retarded Green function can be obtained 
holographically. 

To proceed, we would like to start with the Dirac equation ( p.7|) in a more general 
static background, i.e., 

ds' = -guir)dt^ + grrir)dr^ + 9xAr)[idxY + (dx^], Aa = A(r)(rft),. (3.1) 

Then we choose the orthogonal normal vector bases as follows 

I , d , I , d , 1 



from which the non-vanishing components of spin connection can be obtained as 



(3.2) 



(l^tr)a — ~{<^rt)a 



drx/gti 



{dt)a, {uJi, 



[UJr 



dryMr 



Y 9rr Y 9rr 

Next let tjj = {—h)~iip, then the bulk Dirac equation can be expressed as 

ipT'^^drAt 



-{dx')a. (3.3) 



/9t 



'9x 



'2y/9u9r 



-!f = 0. 



(3.4) 



By the rotation symmetry in the spatial directions, without loss of generality, we shall 



let (fi 



-iujt+ikx 



LP, thus the Dirac equation reduces to 



/9x 



ipV^drA 
^JgTt 



t \ ~ 



)ip + {-iur + ikV^)^ = 



(3.5) 
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where 



u = ^-^{uj + qAt). 



(3.6) 



Set = ( , then with our representation of Gamma matrices, the equation of 
motion can be further simphfied as 

V 9xx 



[dr + m^/g^.a'i)(pi = io2U + [(-1) k 



^1 



(3.7) 



with I = 1,2. Moreover, by y?/ = ), the above equation of motion gives rise to 
the following flow equation, i.e.. 



^9xx „ ^ . ^ r Py/OxxdrAt i \7; i , r , P\^9xxdrAt i -,^2 /o o\ 

= -2mv^6+[M — h(-l) k] + [u-\ — (-1) k]^j (3.8) 



/9r 



where ^7 = |j- Now plug the particular background ( p.2|) into the above equations, we 
end up with 



{r'^/fdr + rma^)(pi = + qgpQil - -)] + [(-l)'A; 



] (3.9) 



for the equation of motion, and 

r^^ffdr^i = -2mrii + + (-l)^fc] + K - (-1)'^]^/ 
for the flow equation, where 

v± = ^[uJ + qgpQil - - ± . 

V/ r r 



(3.10) 



(3.11) 



Whence near the boundary, namely when r goes to the infinity, ipi behave in the 
following way, i.e.. 



(3.12) 



The ratio Gj = ^ can be fixed by imposing the in-falling boundary condition for at 



the horizon, where ipi behave as 



(3.13) 
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with = f Alternatively, this ratio can also be obtained in a more convenient 
way as 

Gi = lim r^""^!, (3.14) 

r— >oo 

by solving the flow equation (|3.10 ) with the boundary condition at the horizon 

= t. (3.15) 



Then by the recipe of AdS/CFT, it follows from the variation of on-shell action ( p. 18 



that at the non-relativist ic fixed point the retarded fermionic Green correlator can be 
extracted from the following relation, i.e., 

with {Di, D2) = + d2, di - 62) and (Ci, C2) = 75(^1 + C2, C2 - ci). To be more 

precise, from such a relation, we have 

7^ + /3 = l, ^-/3 = 1, 
^ + r]Gi = -l, 7-7^^2 = 1, (3.17) 

which gives Or in terms of Gj as^ 

Gr = ^^+2^^ . (3.18) 



^1—02 

G1+G2 Gi+G 



Consequently, det^Gn) = —1, and the trace along with the eigenvalues of Gr can be 
worked out as 



G1G2 — l±>yi + Gl + Gl + G\G\ fn ^c^\ 

a[TG2 ^^-'^^ 



IrGn = A, + A_ = ^^i^. (3.20) 

Lri + Lt2 

It follows from the flow equation ( |3.1U| ) that Gi and G2 are related to each other as 

G2{uj,k) =Gi{u,-k). (3.21) 



"'it is noteworthy that Gi are simply the components of diagonal retarded Green function in the 
standard quantization. An intriguing question is how to understand (3.18) directly from the RG flow 
induced by the double trace deformation. 
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Figure 1: The 3d and density plots of spectral function for the case of q = 1 and p = 0. 

Therefore both of the trace and eigenvalues of our retarded Green function are invariant 
under the transformation k — )■ —k as it should be, guaranteed by the rotation symmetry 
mentioned above. In addition, by the flow equation, we also have 

Giioo, k- gp) = -G;{-uj, k; -gp), (3.22) 

or equivalently 

Gi(w, k; q,p) = -GK-u, k] -q, -p). (3.23) 

To make our life easier, in what follows, we will set gp = ^ and work exclusively with 
the case of m = 0, where the flow equation further implies 

G,iu;,k;p) = --—^ ^. (3.24) 

G2{uj,k; -p) 

So it is essentially enough to restrict ourselves to non-negative fc, q and p. 
3.2 Numerical results 

In the following numerical calculations, we shall focus exclusively on the probe fermion 
in the zero temperature soup, which can be achieved by setting Q = \/3- Then the 
chemical potential is given by = v^g. 

We start by fixing q = 1 but varying p. The corresponding numerical results are 
plotted in Figures 0, |] , |^, and H, where the 3d and density plots are drawn for the 
spectral function on the left and right respectively. First, albeit with the different choice 



Figure 2: The 3d and density plots of spectral function for the case of g = 1 and p = 2. 




Figure 3: The 3d and density plots of spectral function for the case of q = 1 and p = A. 



of Gamma matrices the result presented in Figure |l] is consistent with that obtained in 
ITf , where instead of the spectral function, ImA+ and ImA_ are plotted individually^. 



^As pointed out to us by Robert Leigh, ImA+ and ImA- can also be measured by the so called spin 
polarized ARPES. The reason why we work only with the spectral function lies in the fact that unlike 
the trace of retarded Green function the expression for two eigenvalues generically involves the square 
root, as shown in ( 3.19| ). Such a square root makes it less easy to separate these two eigenvalues by our 
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n 11 
Figure 4: The 3d and density plots of spectral function for the case of q = 1 and p = 8. 



q=l 


p=0 


p=2 


p=4 


p=8 


kp 


No 


1.08188 


2.63542 


4.87894 



Table 1: The Fermi momentum for the case of g = 1. 



This can be regarded as sort of consistency check on our numerics. Second, with 
the finite chemical potential the infinite fiat band get mildly dispersed. In the large 
momentum limit, the band becomes asymptotically flat. This is reasonable as at large 
momenta the corresponding modes sit outside of the light cone and can not decay. 
Furthermore, as one goes to the large momentum limit, the fiat band is always shifted 
to = — v^, independent of the specific value of p, which arises because the frequency 
is measured with respect to the chemical potential. But the larger becomes the value of 
p, the larger becomes the momentum at which the peak of fiat band tends to be sharp, 
which is somehow related to something more interesting occurring in the region of small 
momenta. The infinite band, which is destroyed or depleted from zero momentum 
up to the momentum of the order of the finite chemical potential by colliding our 
probe fermion with the relativistic soup, is gradually recovered to extend down to zero 
momentum as one pushes the bulk dipole coupling up. In particular, at the momenta 
lower than the value of p the flat band disperses the other way than it does at large 
momenta. This implies that it costs more energy to excite small momentum modes 



numerical calculations. The only exception is the case of to = and p — 0, where the two eigenvalues 
can be explicitly massaged as A+ = §7TT ^'^'^ ^ \^%\ using (3.24). 
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Figure 5: The Fermi momentum can be fixed to the 5th digit as kp = 2.63542 by taking 
Lu = —0.00001 for the case of q = 1 and p = 4. 



than intermediate momentum modes, which is also consistent with the phenomenon 
that the peak of flat band becomes sharper and sharper at small momenta as p is 
increased. Moreover, as shown in Figure ^, some other gap states, which are nearly flat 
finite bands, are generated at small momenta when the value of p is large enough. 

On the other hand, as p increases, the Fermi surface starts to show up at some 
point with the Fermi momentum increasing. To be more precise, as demonstrated in 
Figure ||, we can identify the location of Fermi surface to the 5th digit using the fact 
that the location of peak approaches the Fermi surface kp in the limit w — )■ 0. The 
corresponding results are hsted in Table It is noteworthy that such a situation 
is different from what is happening to the relativistic fixed points, where the Fermi 
surface disappears and instead the gap forms when the bulk dipole coupling is large 
enough |jl2| . 



Now let us move on to the case of g = 2. As shown in Figure |, 0, H, and ^ some 
new features come in, namely, the original Fermi surface disappears, accompanied by 
a gap opening up at negative u as one increases the dipole coupling p from zero. But 



^Note that in Figure ^ the apparent disappearance of the peak around a; = is sort of numerical 
artifact. It arises because the peak becomes sharper as one approaches to = 0. 
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Figure 6: The 3d and density plots of spectral function for the case of q = 2 and p = 0, 
where the Fermi surface shows up with the Fermi momentum kp = 1.53521. 




Figure 7: The 3d and density plots of spectral function for the case of q = 2 and p = 2, 
where the Fermi surface disappears and a flat gap sits at w = —0.85. 



as p is increased further, the gap gets bent gradually towards oo = 0. Finally the gap 
closes up, followed by the sprouting of a new Fermi surface. This suggests that there is 
some kind of competition between the charge q and the bulk dipole coupling p. Each 
of them prefers to create their own Fermi surfaces. The gap exists at the parameter 
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Figure 8: The 3d and density plots of spectral function for the case of g = 2 and p = 4, 
where the gap bends towards a; = with the gap energy u = —0.56 for the zero momentum 
mode. 




Figure 9: The 3d and density plots of spectral function for the case of q = 2 and p = 8, 
where the bent gap sprouts a new Fermi surface with the Fermi momentum kp = 1.20996. 

place where they counterbalance each other. Such a pattern can be checked to persist 
for other larger charge cases. Finally, we would like to emphasize that such a rebirth of 
Fermi surface from the generated gap like a phoenix is never occurring at the relativistic 
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fixed points, where instead the gap is generically enlarged as the bulk dipole coupling 
is increased ||T3[|. 

4. Conclusions 

We have worked with one of the recently discovered non-relativistic fermionic fixed 
points and investigated how the corresponding spectral function is modified by the 
bulk dipole term numerically. As a result, although the infinite fiat band is robust 
against the bulk dipole coupling as well as chemical potential, the bulk dipole coupling 
modifies the fiat band in particular at the momenta lower than the value of bulk dipole 
coupling, bending the band towards more negative frequencies. More remarkably, it is 
shown that a new Fermi surface emerges from the gradually disappearing gap as the 
bulk dipole coupling starts to outweigh the charge parameter, which is totally different 
from what is happening at the relativistic fixed points. It is definitely worthwhile to 
identify such features by scanning the realistic condensed matter systems and see what 
kind of role the bulk dipole coupling mimics there. 

We conclude with various directions worthy of further investigation in the near 
future. First, although for the cases of m 7^ the retarded Green function is expected 
to exhibit the qualitatively similar behaviors as the massless story we have told above, 
it is interesting to investigate how some specific features of retarded Green function 
like the location of Fermi momentum depend quantitatively on the mass m. Second, it 
is intriguing to heat up our system to a finite temperature by considering the charged 
AdS black hole away from the extremality to see which features are smoothed out by 
such a finite temperature while which features can persist against the temperature. 
Last but not least, as is well known, the low energy behaviors around the Fermi surface 
are controlled by the AdS2 x region near the horizon for the relativistic fixed points. 
On the other hand, only in how to massage the boundary data does the difference lie 
between the non-relativistic and relativistic fixed points. The bulk dynamics is the 
same for all of these fixed points. So it is expected that at the non-relativistic fixed 
points the low energy behaviors around the Fermi surface are also determined by the 
emergent IR CFT associated with the AdS2 x region somehow. In particular, it is 
highly desirable to have an analytic understanding of the low energy behaviors around 
the Fermi surface at the non-relativistic fixed points by performing the rigmarole of 
matching calculation. We expect to explore these issues elsewhere. 
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